Contraction Maps and Equivalent 
Linearization* 

By J. M. HOLTZMAN 

(Manuscript received July 26, 1967) 

This study is primarily concerned with the question: If the method of 
equivalent linearization indicates the existence of a periodic solution, is 
there actually a periodic solution near the approximation of equivalent 
linearization? To answer this question, ice use a modification of the con- 
traction mapping fixed point theorem. We discuss applications to differential 
equations and difference-differential equations (with forcing functions). 
Also, we show that our use of contraction maps is not applicable (without 
modification) to autonomous systems because the mapping evaluated in 
the neighborhood of a periodic solution to an autonomous system is not a 
contraction in a space of periodic functions. 

I. INTRODUCTION 

The method of equivalent linearization is a most valuable tech- 
nique to investigate nonlinear phenomena, particularly nonlinear 
oscillations. It has its roots in the method of Krylov and Bogoliubov 
and is related to (or equivalent to, depending on the specific defini- 
tions) the method of harmonic balance, Galcrkin's method, and the 
describing function method used by control engineers. The purpose 
of the present study is to develop a new technique for investigating 
the method of equivalent linearization. 

We shall be primarily concerned with the following question: If 
the method of equivalent linearization indicates the existence of a 
periodic solution ,r„. is there actually a periodic solution near .t ? To 
answer this question we first introduce a convenient modification of 
the contraction mapping fixed point theorem which is actually more 
general than just applicable to the question posed above. f We apply 

* Taken from ;i dissertation submitted to the Faculty of I lie Polytechnic In- 
stitute of Brooklyn in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy (System Science). 1967. 

t Appendix A contains some reading suggestions for engineers interested in 
this work but who arc not familiar with the mathematics used. 
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our approach to systems described by nonautonomous differential 
equations. Then we show that there is no essential difficulty in also 
handling difference-differential equations. 

We shall try to clearly indicate what our method can and cannot do. 
The discussion of autonomous systems is particularly important in 
this regard. The relation of the present study to previous work is dis- 
cussed in Section VIII. 

II. THE METHOD OF EQUIVALENT LINEARIZATION* 

Consider the following vector differential equation 

*(0 = Kx(t), i) = A(t)x(t) + n(i, x(t)), (1) 

where 

f(x, t) = f(x, t + T) (2) 

for all {x, t) of interest. This, of course, includes the case of j(x, t) 
independent of t, i.e., 

Kx, t) = f(x). (3) 

We shall be concerned with the situation that permits an equivalent 
representation of (1) : 

x = LN(x), (4) 

where x now represents a vector function, L is a linear operator, and 
A^ is a nonlinear operator (these terms will be made more precise 
later). If it is assumed that LN{x) has the following Fourier series, 

LN(x)(l) ~ | + Z (a, cos k y t + h sin k ^ ') ■ (5) 

Then we define LN(x) as follows: 

LN(.r)(t) = a, cos -= t + b, sin ^r t . (G) 

That is, L extracts the fundamental component of the Fourier series. 
The method of equivalent linearization seeks a solution of the equation 

x = LN(x). (7) 

This study will be primarily concerned with the following problem: 



* See Minorsky, 1 p. 350, for a discussion of the relationship of the method of 
equivalent linearization to the method of Krylov and Bogoliubov. 
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Given an x satisfying (7), is there an x* satisfying (4) and if there is, 
how are x , and x* related? 

Note that (4) is a functional relation more general than (1) and 
our method will be correspondingly applicable to a more general 
problem. 

The above discussion is now related to the method of describing 
functions* as commonly used by control engineers. They are con- 
cerned with the feedback loop shown in Fig. 1. The linear operator L 
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Fig. 1 — Feedback Loop. 

is represented in this case by a transfer function H(u>) (see Kaplan 4 
for a definition and discussion of transfer functions) and the nonlinear 
operator N is represented by a nonlinear function, 

y = n(fii - x), 

where ej is an input function. The engineer replaces the nonlinear 
function n by its describing function which is defined loosely as the 
complex ratio of the fundamental component of the output to a 
sinusoidal input. That is, iff 



n(A sin wt) = ^2 a k sin hat + b k cos kut [u = — 
then the describing function of n is 



Vai + 6? /tan" 1 ^ 



Note that while the describing function may be dependent on both A 
and a), it is still a relatively simple matter to replace n by its describ- 

* For further discussion of the use of describing functions by engineers see 
e.g., Truxal 8 or Graham and McRuer.' 1 They give further references and historical 
background. The describing function method is associated with the names of 
Tustin, Goldfarb, Oppelt, Kochenberger. Dutilh, and Nichols and Krcezer. Also 
see Minorsky, 1 Chap. 17 for a discussion of the work of Theodorchik and Rla- 
quiere. The work of E. C. Johnson is discussed in Ref. 2. 

t The constant term is assumed zero. 
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ing function, then to consider it as a "linear" (or quasi-linear) op- 
erator and then use standard techniques for linear systems. Of course, 
such a procedure should be mathematically justified and, in fact, that 
is the purpose of this study. 

Before we embark on our investigation, it is well to review the argu- 
ments used by engineers in their justification of the method. These 
arguments seem to be plausible and they are suggestive of what may 
be expected of a more rigorous investigation. If it is assumed that the 
combination of n and H{w) operating on a sinusoidal function is pri- 
marily fundamental (i.e., the harmonics are "small" compared to the 
fundamental) then it would be expected that the describing function 
method might not be too inaccurate. The harmonics will be small if 
one or both of the following are satisfied: 

(i) the nonlinearity n is "not too nonlinear" 

(ii) the transfer function H{oi) is low-pass, i.e., it attenuates har- 
monics much more than the fundamental. (It is assumed that no sub- 
harmonics arise). 

We shall use Duffing's equation, 

ij + ay + by 9 = ! cos ul, 

as a running example to illustrate the methods discussed. We show 
here how this differential equation corresponds to a feedback control 
problem and then make no further explicit reference to feedback sys- 
tems. The appropriate feedback system is shown in Fig. 2. 

The next section contains an approach to a problem much more 
general than the problem of equivalent linearization posed in this 
section. The remainder of the study will be primarily devoted to 
adapting the more general approach to the specific problem of equi- 
valent linearization. 

It may be noted that we are not getting more abstract in the next 
section just for the sake of abstraction. It should be clear to the reader 
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Fig. 2 — Feedback equivalent of Duffing's equation. 
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that the method of equivalent linearization leads to an integro-differ- 

cntial equation rather than an ordinary differential equation because 
Fourier coefficients are determined by integration (this is also pointed 
out in Bass 8 , p. 898). We then cannot expect the theory of ordinary 
differential equations to answer our questions and we are led quite 
naturally to considering more general equations. In particular, the 
theory of operator equations in a Banach space is shown to provide 
the tools appropriate to the task. As an added bonus for the abstrac- 
tion, wc develop an approach which is applicable to problems unre- 
lated to equivalent linearization. 

III. THE USE OF THE CONTRACTION MAPPING FIXED POINT THEOREM WITH 
DERIVATIVES IN A BANACH SPACE 

Let X be a complete metric space (with metric d) containing the 
the closed set n and let P map Q into itself. P is a contraction map- 
ping if 

d(P(x), P(x')) ^ ad(x, x') (.r, x' t 2) (8) 

with a < 1. The contraction mapping theorem t states that if P is 
a contraction mapping then there is a unique x* e Q such that x* = P(x*), 
i.e., x* is a fixed point of the operation P. x* is the limit of a sequence 
j.r„i whore 

•r n+1 = P(x„) (9) 

and Xo is any element of Q. Furthermore, 

<*(*. , x.) £ ***** = ''^• x " ) ,, = 1,2,.... (10) 

1 — a 1 — a 

In order to use the contraction mapping fixed point theorem it has 
to be shown that some neighborhood of .r is mapped into itself and 
that in this neighborhood the operation is contracting. Our approach 
will be simultaneously to determine a set containing .t which is mapped 
into itself along with the contraction constant a for the operation on 
that set. This is possible because of relationship (10). The use of opera- 
tor derivatives will be seen to be convenient. The method will result 
in a relation in a for which it is desired to find solutions with a e [0, 1). 
The following is proven in Kantorovich and Akilov (p. 601). 



t Sp(! Kiintorovicli mid Akilov." ]«. (i'27. 
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If X is a Banach space and P maps a convex closed subset ft of X 
into itself and if P has a derivative* at every point of ft, then 

sup || P'(x) || = a < 1 (11) 

scD 

implies that P is a contraction on ft (and thus, there is a unique fixed 
point of P in ft). 

The object then is to find a neighborhood of * mapped into itself 
and in which the norm of the derivative is less than one. The follow- 
ing simple theorem is a help in this direction. 

Theorem: Let B be a Banach space. F maps B into itself and x t B. 
It is assumed that 

(i) F has a derivative at all x e B 
(ii) There is a nondecreasing junction g such that if x t B, then 

\\F'(x) || ^ g(\\x-x Q ||) 

(Hi) There is an a c [0, 1) such that 

k 



9{—J * «> 

where 

k £ \\F(x Q ) - x || 

Then there is a unique x* e ft such that 

x* = F(**), 

where 

k 



ft = \x:xzB, || 2 - x || £ j _ 

Proof: We will show that || F'(x) || ^ a for all x e ft and that F maps ft 
into itself and thus, there is a unique fixed point in ft. If x c ft, we have 
from (ii) 

\\F'(x) || Srfll*-*. ID 

I: 



\1 — or 



< 



*See Kantorovich and Akilov,' chap. XVII, for a general discussion of dif- 
ferentiation in Banach spaces. For convenience, Appendix B of this study repeats 
the definitions. 



EQUIVALENT LINEARIZATION 2411 

F maps 12 into itself because if x t 12, then 

II F{x) - x || £ || F(.v) - F(x ) ll + H F(x )- x || 

£a a 1 1 x — x 1 1 ;+ k 

^ a y^— + k 
1 — a 

k 



1 - a 
We further discuss the use of the contraction mapping theorem in 
Ref. 7. This modification of the contraction map theorem is less gen- 
eral than some other modifications but is simpler to apply when ap- 
plicable. 

IV. THE CHOICE OF BANACH SPACE 

In order to use the result of the previous section, an appropriate 
Banach space must be chosen. For most of our investigation it will 
be found convenient to use the space of continuous periodic functions. 
Another space worth considering is the space of periodic functions 
square-integrable over a period. Before discussing the desirable char- 
acteristics of this space, a restrictive factor will be mentioned. The 
nonlinear operator of interest y — N(x) is often defined by the ordi- 
nary function (satisfying the Caratheodory condition) 

7/(0 =n(t,x(t)). (12) 

A necessary condition that this operation map L 2 (0, T) into L 2 (0, T) 
is (see Krasnosel'skii, 8 p. 27) that for some b > and some a(t) c L 2 (0, T) 

I n{t, u) | ^ a(t) + b | u \ t e:[0, T]. (13) 

It is thus seen that the allowable nonlinearities are quite restricted. 
This is, in fact, the reason the present investigation will be carried out 
in a space of continuous functions where the requirement that a 
function map a continuous function into a continuous function is 
much more convenient. It should be noted, however, that in some 
cases one may focus attention on some subset of the Banach space 
and less restrictive requirements on the nonlinearities might be im- 
posed. Also, for many control engineering problems the nonlinearities 
are Lipschitzian and the problems can be attacked in Lg. 

The attractive feature of L 2 is that Fourier series results can be 
fully utilized fin particular, Parseval's relation). More generally, L 2 



2412 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1967 

is a separable* Hilbert space with many useful properties and the 
trigonometric functions arc a complete orthogonal system in L 2 . The 
norms can often be conveniently evaluated in terms of quantities as- 
sociated with the "transfer function" or "frequency response". For 
example, L may be defined by the set of complex numbers {. . . , L 2 , 
L_j, L , Li, L 2 , . . .} (i.e., the transfer function evaluated at the fun- 
damental and harmonic frequencies) .f A simple sufficient condition for 
L to map L2 into itself is that 

sup | L n | < oo . 

The evaluation of || LN{x Q ) - .r„ || may be done as follows: 
|| LN(xJ - ;r„ || = || LN(x ) - LN(x ) \\ 
£ || L - L ||-|l N(x n ) H 
|| L || = sup ] L„ |. 

The last relationship is proven in Appendix B of Sandberg. 10 

Despite the above mentioned attractive features of L 2 , we chose 
to work in the space of continuous functions primarily because of the 
first-mentioned restriction placed on the nonlinearities in L>. Also, 
the sup norm (uniform norm) in the space of continuous functions 
seems more appropriate in error analysis (the error between an ap- 
proximation and an exact solution) than does the L 2 norm. The sup 
norm provides a bound on the magnitude of the error while the L 2 
norm gives the integral of the square of the error. 

Section V will give the details of working in the space of continuous 
periodic functions. First (in Section 5.1) an integral equation equi- 
valent to the differential equation of interest will be derived. Then 
in Section 5.2, the derivatives will be determined and finally in Sec- 
tion 5.3, the quantity || F(x ) - x || will be evaluated. Application of 
the results will then be seen to be rather straightforward. 

V. APPLICATION TO DIFFERENTIAL EQUATIONS 

5.1 The Equivalent Integral Equation 

Halanay 11 shows how to convert the quasi-linear differential equa- 
tion 



*It is, of course, assumed that (ho measure is Lebesgue. Then U is separable; 
see Kolmogorov and Fomin," Vol. II. p. 88. 

t NJ confusion should arise because of the double use here <>l the symbol L : , 
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dx(t) 



(It 



= A(f)x(t) + n(t, .c(0) (14) 



into an integral equation which is convenient for examination of pe- 
riodic solutions. First, consider 

^ = A(t).v{l) + /(/), (15) 

where .4(0 and fit) are both continuous and periodic of period T. 
The following theorem is proven in Halanay, 11 p. 223. 

Theorem: A necessary and sufficient condition i?i order that, jor any 
periodic junction f{t) of period T, system (15) admits periodic solu- 
tions of period T is that the corresponding homogeneous system 

<hl ( ^ - A{t)y{t) (16) 

does not admit a non-trivial periodic solution of period T. 

It is to be noted that if Y[t) is the principal fundamental solution 
matrix for (16) then the existence of the inverse of [/ — Y(T)]* is 
equivalent to the non-existence of a non-trivial periodic solution (of 
period T) to (16). Then the following proposition is proved in 
Halanay," p. 225. 

Proposition: If \I — Y (7')| ' exists, the unique periodic solution of 
the system (15) can be put in the form 

•'•(/) = / 0(1, s)f(s) ds, (17) 

where 

^-j^-wrr-^), os.*fsr (18) 

{Y(t + T)[l - Y(T)]- 1 )-'(s), ^ t < 8 ^ T. 

Since this reformulation into an integral equation is quite important, 
a sketch of the proof given in Halanay 11 will be given here. Solution 
of (15) is 

,•(/) = Y(t)x(0) + f nt)Y-\ S )1(s) ds. (19) 



* / is the identity matrix. 
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For periodicity, 

x(T) = x(0) = Y(T)x(0) + f Y(T)Y-\s)f(s) ds (20) 

Jo 



or 



[/ - Y(T)]x(0) = f Y{T)Y-\s)j{s) ds. (21) 

Jo 

Since [I - Y{T)] is assumed invertible, we can solve for x(0) and get 
x(t) = Y(t)[I - F(r)]- 1 f Y{T)Y-\s)j{s)ds 

Jo 

+ f Y(t)Y-\s)Ks)ds. (22) 
Jo 

The form of G(t, s) given in the statement of the proposition results 

from algebraic manipulation of (22). 

Now suppose that 

^ = A(t)x(t) + n(t, x(t)) (23) 

at 

with A{t) and n(t, u) both being periodic of period T and [/ - Y(T)] 
invertible. It is assumed that n{t, x(t)) is continuous if x{t) is con- 
tinuous. Then, from the previous discussion, if we can find a continu- 
ous periodic x of period T satisfying 

x(l) = f G(t, s)n(s, x(s)) ds, (24)* 

Jo 



we have a periodic solution of (23). (It is easily shown that such an 
x satisfies (23) ; see Halanay, 11 p. 237) . 

The problem is thus reduced to finding a continuous solution of a 
nonlinear integral equation. We need only consider the interval [0, T] 
because G(t, s) was constructed so that .r(0) = x{T). 

5.2 Computation of Derivatives 

It is shown here how to evaluate the Frechet derivative! of the 



* The nonlinear integral operation represented by the right hand of (24) is of 
the form sometimes referred to as a Hammerstein operator which is a special case of 
Uryson's operator denned by /J" K(t, s, x(s)) ds (see Krasnosel'skii, 8 pp. 32, 46). 

t The Frechet derivative is actually more than what ia required. The Gateaux 
derivative (which does not require uniform convergence) would suffice for much 
of what follows. However, since the convergence is indeed uniform in most cases 
of interest and since the uniformity is easy to demonstrate, we shall derive the 
Frechet derivative. Furthermore, Frechet derivatives are needed in Section VI. 
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mapping y = F(x) defined by 



//(/) = f G(t, S )n(8, x(s)) ds i t [0, T]. 



(25) 



This operation is assumed to map into itself the Banaeh space of 
real-valued n-veetors continuous on [0, T] with norm 



II* 



= max 



max 

tt\0,T] 



x t {t) 



(26) 



To determine the derivative of the operation y = F(x) it is con- 
venient to express it as y = LN(x), where N(x) is defined by the 
nonlinear function n(t, x(t)) and L is the linear integral operator. 
Then F'(x ) = LN'{xq) (see Kantorovich and Akilov, 6 p. 659). Thus, 
consider the mapping y = N(x) defined by 

y(t) = n(t, x(t)) 

fn,(/, x(t)j] 

: • (27) 

It is assumed that n(t, x(t)) is continuous whenever x(t) is continu- 
ous on [0, T]. For simplicity, the derivative will be determined for 
the case of 



(*,*(*)) =0 i = 1,2, ••• ,n - 1 



(28) 



n n (t,x(t)) = p(t)h(x 1 {t)) +r(t), 
where p(t) and r(t) are continuous functions of t with period T and 
h(u) is a twice continuously differentiate function of u. This special 
case which covers our examples may arise, for example, when the 
matrix differential equation is actually derived from a scalar differen- 
tial equation. The more general case offers no other difficulties than 
much more complicated notation (e.g., one must deal with matrices 
of partial derivatives). 

The derivative operation * = N'(x )x is defined by 



*(0 = 











pWixoMxi®. 



(29) 
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where 

X (t) = (.1-0,(0, ■ • • , X 0n (t)) T , •'•(/) - (*i(«), ■ ' ' ! *»(*)) r , 

(superscript T denotes transpose) and 

h'{xM) = -^- 

To prove this, it must be shown that 

lim N(x ° + ^ - y w = pw., 

that is, 



(30) 



(31) 



lim max 

,1-0 ItlO.Tl 



p(t)h(x ni (l) + M.r,(0) - p(t)Hx ni (t)) 

- p(t)h'(xM)xi(t) 



- 0, (32) 



with that convergence being uniform with respect to all .r with || x || — 
1. Since 



p(0^oi(0 + gaga - p(t)h(x 0l (t)) _ p{t)h , Mt))Xi{t) 



= I p(0 



h(x 01 (t)n^(t)h'(x 01 (t)) + ^^ fc'WO+*Ol*i(0)--*(*»W) 



= I Pit) 



- h'(x 0l (t))xi(t) 



■ | .r,(/) | 2 | h"(x 0l (t) + Hfin*M 



< 0(0 < 1 



max | p{t) | max | h"{z) | (for | y. \ S 1 and || * || - 1), 

/r(O.T| :tZ 

(33) 



where 



Z = |z:z = w + »;w = x in (t),tt[0,T];\v \ £ 1} (34) 

the uniform convergence relationship is seen to be satisfied. 
To summarize the result, the first derivative operation y = F'(x )x 
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is defined by 



,v(o = r G(t, s ) 



ds. (35) 



jK«)V(a.l(«))*i(«) 

Loosely speaking, the derivatives of the integral operators are ob- 
tained by differentiating under the integral. 

5.3 Evaluation of || LN(x„) — x || 

A bound on || LN <.r„) - x || is given here. The relationship 

\\LN(x ) -.fo|| ^ \\LN(xo) || + ||so II ( 36 ) 

is too gross an estimate. The evaluation is simplified if the following 
relation is used : 

LN(x ) - x = (L - L)N(x ). (37)* 

Recall that the operation L suppresses all frequency terms except 
the fundamental. 

Consider the same system as in Section 5.2 and assume that 

p(t)h(x t (t)) + r{t) = f] a K cos fc|i+ b K sin k y t. (38) 

and also, for simplicity, that .4 is a constant matrix. Then, 
\\LN(x ) - x || 

C T X ( 9 7T ') T \ 

= max max / G in (t, s) £ [&k cos k 7F s + b * sin k V s ) ds 
^ E (I <*k I + I h K |) max max f | G in (t, s)\ds. (39) 

fe_ 2 1-1, ••*." ttlO.Tl •'O 

5.4 Example 

Consider Duffing's equationf 

y 4- ay + &y a = / cos cof (a > 0). (40) 

Equivalent linearization indicates that 

y = A cos o>t (41) 



*From this expression, it is seen that \\LN(x a ) — x ||may be regarded as a 
quantitative measure of characteristics (i) and 07) mentioned in section II. 

t DuHing's equation is discussed in great detail in Stoker.'- Also sec Graham and 
McRuer 1 for a treatment of DuHing's equation as a feedback control problem. 
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with 

IbA 3 + (o - co'V - / - 
is an approximate solution of the equation. 
Letting 



1007 






X 

x 2 



.}-. = y 

.r._, = Xi 
the corresponding vector differential equation is 



x = 



x + 



1 

l-a OJ 
The fundamental solution matrix for 

1 

-a 




/ cos cot — bxy 



y = 



IS 



Y(t) = 



COS Yd t 



in Va t 



sin 



G{t, s) is given by 
G(t, s) 




Va 
— Va sin Va £ cos Va £ 



(42) 

(43) 

(44) 
(45) 

(46) 

(47) 

(48) 



sin Va (--f+s) ^ cos Va ^2 _i + s 

- Va cos Va («—*+«] sin Va ^- —t+s) 



—sin Va (— +£— s 



^ s ^ J ^ T (49) 



= cos Va (— +/— s 



Va 



- Va cos Va (-+<-s) -sin Va (j"+'~ 8 J 
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The integral operation of interest, y = F(x), is defined by 

" ds. (50) 

./ cos cos — bx\{s), 

The approximate solution (obtained from harmonic balance) is x , 
i.e., 



ij(t) = f G(t, s) 

•'o 



■Vo(t) = 



U 02 (OJ 



A cos ut 
—uA sin at. 



(51) 



with co and A being related by (42). 
From Section 5.3 we have that 



LN(x ) — -T 1 [ ^ T max max | G i2 (t, s) 

1-1.2 I ,ac[0,T\ 



bA- 



< TC 



bA : 



There 



C = 



sin 



VnT 



max {1, 1/Vo), 



(52) 



(53) 



.(0 = f G«, •) 

•'o 

-I 



ds 



ds. 



(54) 



The derivative operation z = F'{xq)x is given by (see Section 5.2) 


[— 3fc.roi(s).T,(s). 
~-<r' 12 (/,s)36.ii(s) dfxiisj 
-G 22 (t, s)3bxl(s) oj .r 2 (s) j 

The norm of the derivative operation at an arbitrary point x (not 
necessarily at .r n as above) is evaluated as follows: 

1 1 F'(x) 1 1 ^ max max f \ G i2 {t, s)Sbxl(s) \ ds 

i = 1.2 ItlO.T] -'o 

g 3 | b | CT( max | Xl (t) \f 

itio.n 

^ 3 | 6 j CT( max | .r (11 (/) | + max | x t (t) - x ni (t) |) 2 

UIO.T] (t[o.7'| 



^ 3 | b | 6'7 T (| A 1 + 11* -.(•„ 



(55) 
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The above relation defines the nondecreasing function g such that 



|| F'(x) || ^ p(||* - x, 
to use the theorem of Section II, let 



il = <x: II .r - .r„ || ^ z 



1 - a 



where 



k = CT 



\bA : 



(56) 



(57) 



(58) 



(see (52)). 

If an a e [0, 1) can be found satisfying 



3 I & I CT 



AI + 



CT 


&A 3 |] 


4 


1 


— a 



< 



(59) 



then there is an x* c il such that :c* = F(x*), i.e., Duffing's equation 
has a periodic solution in the neighborhood of the approximation obtained 
by harmonic balance. 

Rather than just solve the cubic relation (59) for particular numerical 
values of a, b, /, and co (which is, of course, the thing to do if one is 
given a particular equation of interest) we shall obtain some general 
results. Consider a, f, and w fixed and / ^ 0. Sincet for any « t |0, 1) 



lim 3 | b | CT 
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(60) 



it is seen that for b sufficiently small, there will be an a t [0, 1) satis- 
fying (59) (and thus a periodic solution neighboring the approxima- 
tion). Note that while the result has been stated as an asymptotic 
result, it is possible to determine quantitatively what is meant by 
"sufficiently small". This is in contrast to most asymptotic analyses 
based on "small nonlinearities". 

5.5 Special Cases 

In many cases, it is not necessary to convert the differential equa- 
tion into a vector integral equation. For example, let the system be 



f See Appendix C lor details, 
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described by the equation 

a () .r (w, (0 4- ay-^il) + • • • 4- a n x(t) = n(t, x(l)), (61) 

where the a„ n,„ arc constants, n(t , x) is continuous in t when 

.r is continuous in t, and 

n{t + T, a) = n(t, u). (62) 

A periodic solution to the differential equation will satisfy the follow- 
ing integral equation 



.,•(/) = / I !',.(/ - u)n(u, x(u)) du, 

Jo 



(63) 



where W T (t — it) is the appropriate convolution kernel (see Kaplan, 4 
chap. 4 for details). All the manipulations of previous sections will 
be somewhat simplified as a result of not having to deal with matrices. 
In particular, the example of Section 5.4 could be repeated with some 
simplification. We omit the details because they exactly parallel the 
previous case. We felt it would be more useful to work out the details 
of the more complicated case. As the order of the differential equation 
increases it clearly becomes more advantageous to avoid the use of 
matrices. 

It may be noted that there is a finite Fourier transform Y(iruo) 
(again see Kaplan, 4 chap. 4) associated with the differential equa- 
tion (61) 

Y(inu) = ,. * z (64) 

«„(imw) + • • • + a,„ 

This Y(p), considered as a function of a complex variable p, evidently 
can only have poles and cannot have finite zeros. In many electrical 
engineering applications (e.g., control systems, networks) the relevant 
transfer function has both poles and zeroes. In these cases, we would 
start with the transfer function, rather than a differential equation of 
the form (61), find the corresponding convolution integral and then 
apply our method. For other applications it must, of course, be verified 
that the appropriate conditions are satisfied. 

5.6 Autonomous Systems] 

The describing function method has been used by control engineers 
primarily for the prediction of self-oscillations (i.e., with no forcing 



t Since the actual oscillation of an autonomous system may have a different 
period than that of the approximation, it is usually convenient to normalize the 
time variable and have the period be a parameter. 
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functions). It would seem at first glance that our approach should be 
appropriate for analysis of this problem. Suppose the describing 
function method indicates that there exists a non-trivial periodic 
solution, .To, to the operator equation 

.r = LN(x ). (65) 

Usually, JV(0) = 0, so that it is of interest to investigate whether 
there is any non-trivial solution to the exact equation near x . If our 
method is successful, then we can guarantee that n does not contain 
the trivial solution (x = 0) if 

l l*W-*l l<||»|| (66) 

I — a 
since the fixed point x* satisfies 

i hc . _ ... ii < u F y -"ii . (67) 

Unfortunately, an attempt to use the approach in the autonomous 
case will be unsuccessful. The reason for the failure of our approach 
is due to the nature of the fixed point: the mapping is not a contrac- 
tion in a neighborhood of the fixed point.f The discussion below will 
clarify this point. 

Assume that the differential equation of interest is 

x(t) = Ax(t) + n(x{t)) (68) 

with A a constant real valued matrix and n(x) is a real-valued func- 
tion having continuous partial derivatives with respect to all of the 
elements of the vector x. Suppose that there is a continuous periodic 
x* of period T satisfying (68). Then x* satisfies the equivalent 
integral equation : 

x*{t) = [ G{t, s)n(x*(s)) ds. (69) 

The equation of first variation corresponding to (68) is 

y(t) = Ay(t)+I^y(t), (70) 

where BnJBx* is a matrix with entries dni/dx-j evaluated along the 

t Actually, this should not be surprising since if x(t) is a periodic solution, then 
so is x(t + c). 
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trajectory defined by x*(t). As is easily shown (by differentiation; 
see, e.g., Hochstadt, 18 p. 251) the time derivative of x*(t) satisfies the 
equation of first variation and thus the equivalent integral equation 

x*(t) = J' G{t,s)~i*(8)d8. (71) 

Now consider the derivative operation y = LN'(x*)x associated with 
the integral operation of (69). The derivative operation is defined by 

y(t) = j[* (}{t,s)^ i x{s)ds. (72) 

We have shown above that x* satisfies this last equation or 

x* = LN'(x*)x*. (73) 

Then 

|| x * || ^ || LN'(x*) 11-11** II- ( 74 ) 

Since || x* \\ > 0, we must have that 

\\LN'(x*) || ^ 1 (75) 

or LA' cannot be a contraction in a neighborhood of x*. We leave open, 
however, the possibility that the contraction mapping theorem might 
be applicable in a subspace. 

The above reasoning also shows that there would be difficulty asso- 
ciated with using Newton's method for the problem. To seek a zero 
of P(x), Newton's method uses the following iteration: 

•r n+I = .f„ - [P'(.x n )T\P(x n )) n - 0, 1, ■■• . (76) 

Letting 

P(x) = x - F(x) (F(x) = LN(x)), (77) 

we are led to investigating the inveritibility of I — F'{x) where I is 
the identity operator. Consider the operation 

y = P'(x*)x = (I - F'(x*))x. (78) 

If x = 0, then y = 0. But because F'(x*) is associated with the equa- 
tion of first variation there is also a nonzero x (the time derivative 
of x*(t)) which results in y = 0. There is thus not a unique x satisfy- 
ing y — and P(x*) is not invertible (see Kantorovich and Akilov, 6 
p. 168) . 
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5.7 Properties o] the Fixed Point— Dependence on Parameters; Stability 
The method of this investigation is to obtain a relation for a in terms 

of parameters of the differential equations. If certain conditions are 
satisfied and j£ a < 1, then a represents a contraction constant. 
If a is a contraction constant and, also, a depends continuously on a 
parameter then for sufficiently small changes of that parameter a 
will still be a contraction constant and a periodic solution is still guar- 
anteed. Once again note that "sufficiently small" can be quantitatively 
determined if one wishes to do that. As an illustration, in the Duffing 
equation example, a depends continuously on b. 

If our method indicates the existence of a periodic solution of period 
T then there is no neighboring solution of period T (for a sufficiently 
small neighborhood). This condition does not imply stability of the 
periodic solution. With a perturbation of the initial conditions, stabil- 
ity is concerned with closeness to (asymptotic stability is concerned 
with the eventual approaching of) the original periodic solution. A 
perturbation of the initial conditions of the period solution of period 
T may result in a solution not of period T and thus, not even con- 
sidered in the Banach space used. 

A sufficient condition for the asymptotic stability of a periodic 
solution (of period T) to a nonautonomous system is the asymptotic 
stability of the null solution of the corresponding equation of first 
variation (Hochstadt, 13 p. 251). We are only able to show that the 
equation of first variation may not have a (nontrivial) periodic solu- 
tion of period T if the fixed point is a contraction. The nonexistence 
of a periodic solution to the equation of first variation is (along with 
a continuous differentiability requirement) a sufficient condition for 
the continuous dependence of the periodic solution on a parameter. 
This result is not identical to but is compatible with our initial com- 
ments on the continuity of the contraction constant with respect to a 
parameter. 

5.8 Perturbation Analysis 

A very common approach to nonlinear problems is to solve a linear 
problem ignoring the nonlinearity and then to use a series expansion 
or a perturbation about the linear solution (see, e.g., Hochstadt, 13 Sec- 
tions 6.5, 7.4). As useful as these procedures are, they usually suffer 
from the defect of not providing adequate quantitative information 
about the nonlinear solution, i.e., it may not be possible to determine 
quantitatively what is meant by "sufficient small". Our use of 
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the contraction mapping theorem may prove useful in this regard. 

As an illustration, again consider Duffing's equation and assume 

that Xq was obtained by ignoring the nonlinear term (by 3 ). In this 
case, x = F(xo) is defined by 



(0 = / f*i«(£) s)/ cos ess ds. 



(79) 



Then arguments similar to (but simpler than) those of Section 5.4 
show that if b is "small enough" there is a continuous periodic solu- 
tion to (40) which is "close" to the linear approximation 

... / cos ul ,„_. 

zoi(0 = -'- 2- (80) 

a — co 

Note that "small enough" and "close" may be quantitatively eval- 
uated. 

VI. DIFFERENCE-DIFFERENTIAL EQUATIONS* 

Our use of the contraction map fixed point therom is not limited 
to ordinary differential equations or integral equations. As a further 
example, consider the difference-differential equation represented by 

x = LN(D h x), (81) 

where y = l),,x is defined by 

y(t) = x(t - h). (82) 

If the Banach space B of interest is the space of continuous periodic 
functions, then 

II A, || = 1. (83) 

This follows easily from 

|| D h || =sup \\\D h x \\:xzB, \\x \\ = 1J 

= sup {max | x(t — k) |: x t B, max | x{i) \ = 1} . 

Assume that N maps B into itself. If N is diffcrentiable (i.e., has a 
Frechet derivative) at .r n then M(D h x) has a derivative at x (Kan- 
torovich and Akilov/ 1 p. 658) given by X'{D h x )D„. Then LX(D h ) 
has a derivative at x„ given by LX'(D h x l ,)Dj l . The norm is easily 

* An interesting discussion of (he problems of oscillations in difference-dif- 
ferential equations is given in Chap. 21 of Minorsky. 1 Halanay" contains much 
information (and references) on difference-differential equations. 
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evaluated: 

1 1 LN'{D h x )D h || S 1 1 £ || • 1 1 N'(D h x Q ) \ \ . (85) 

Conceptually, the introduction of the time delay offers no great dif- 
ficulty as compared to the case without the time delay. However, 
it will generally complicate the arithmetic involved in examples, in 
particular, in obtaining the solution to x = LN(x ). This relative 
lack of complication going from differential equations to difference- 
differential equations is not typical. In existence, uniqueness, and 
stability considerations one must consider initial function conditions 
in difference-differential equations while the initial conditions for dif- 
ferential equations are merely at one time (or perhaps boundary con- 
ditions at several times). 

To illustrate the above remarks, consider the following difference- 
differential equation : 

y + ay + hyl = / cos ut (80) 

y h (t) - V(t - h) (87) 

This is Duffing's equation but with the argument of the cubic term 
retarded. The corresponding operator equation is 

x = LN(x) = LiN^Dnx) + F], (88) 

where y = Ni (x) is denned by the cubic nonlinearity, L is the same 
linear operator as in Section V, and 





m = 

It is clear that 



j cos co£. 



(89) 



\\LN'(xo)\\£ || L || || Ni(xo)\\ (90) 

and that the analysis will be completely analogous to that of Section 
V, except that the approximate solution, x , will be different. Note 
that the Banach space is the space of continuous periodic functions, 
not the space of functions continuous on one period. 

To obtain the equivalent linearization approximation let 

y(t) = A cos ut + B sin wt ,»-.. 

= Csin (cot + 9), 
where 
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° " ^^ (92) 

Substitution of this function into (86) yields 

-w 2 C sin (w* + 9) + aC sin («d + 0) + |C 3 6 [cos «fc sin (col + 9) 

— sin coh cos («/ + 0) -f- third harmonics] = f cos «/. (93) 

The approximation is obtained by neglecting the third harmonics and 
equating coefficients of cos cot and sin cot. It is interesting to compare the 
equivalent linearization solution obtained for the difference-differential 
equation with that obtained for the following differential equation 
(Duffing's equation with a damping term) : 

y + ky + ay + by 3 = / cos at. (94) 

Substituting (91) into (94) yields 

~co 2 C sin (cot + 0) + ftwC cos (wi -(- 0) + aC sin («/ + 0) 

+ f&C 3 sin (w/ -I- 6) + third harmonics = / cos cot. (95) 

Comparing (93) and (95) it is seen that, as far as harmonic balance 
is concerned, the effect of the lag is to introduce a damping term 
with damping coefficient k, 

k = -(sinwfc)fC 2 6/co (96) 

(Also, one other term is multiplied by cos coh). 

For some parameter values, the equivalent damping is negative. 
Because of the negative damping, it appears that the periodic solu- 
tion is not asymptotically stable. We say "appears that" rather than 
making a more definite statement for the following reason. While it 
seems plausible that the stability properties of the solution of the 
equation of equivalent linearization should carry over to the actual 
solution, the mathematical proof does not seem so obvious. 

VII. RELATION TO PREVIOUS WORK* 

As mentioned previously, the method of equivalent linearization 
has its roots in the method of Krylov and Boguliubov. For an ac- 

* The literature on equivalent linearization is vast. We shall thus discuss only 
those references which seem most pertinent. Even in those cases, we shall discuss 
only those aspects which are directly related to the present study. The reader 
should consult these references for many other interesting ideas. 
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count of the very important work of Krylov, Boguliubov, and Mitro- 
polsky in this area, see Minorsky. 1 Their work is primarily of the 
asymptotic type, i.e., leading to statements of the form, "for suf- 
ficiently small /x, there exists . . . ." We may view our approach as 
using a fixed point theorem to be able to determine quantitatively 
what is meant by "sufficiently small" for a somewhat different but 
related problem. 

Bass 5 considers the justification of the method of equivalent lineari- 
zation in the autonomous case. In view of our comments on the inap- 
plicability of the contraction map fixed point theorem, it is of interest 
to note that Bass uses a much more sophisticated fixed point theorem. 
Much of his analysis is interesting and important but his final results 
are unfortunately difficult to apply (as Bass himself points out). 

Sandberg 10 considers the operator equation*, 

x = LN(x + j) (97) 

and the equivalent linearization approximation 

x = LN(x + f). (98) 

Sandberg's analysis is carried out in the space of periodic functions 
square integrable over a period. He presents conditions under which 
there exists a unique periodic response to an arbitrary periodic input 
with the same period as well as an upper bound on the mean square 
error in using equivalent linearization. He also gives conditions under 
which sub-harmonics and self-sustained oscillations cannot occur. 
Sandberg's method is to determine conditions that guarantee that 
LNf is a contraction mapping in the whole space. As mentioned pre- 
viously, we do not try to obtain a contraction mapping in the whole 
space but only in a neighborhood of x . We thus free ourselves from 
Lipschitz type requirements. It may be noted that many nonlinearities 
encountered in engineering are non-differentiable and Lipschitzian 
(e.g., piecewise linear functions such as saturation-type nonlineari- 
ties). For these, Sandberg's analysis is applicable while ours is not 
because we have required differentiability. Thus, Sandberg's work 
and ours complement each other in this regard. Also, Sandberg very 
fruitfully uses Fourier transform results in his analysis of feedback 
systems. 

* This is the same notation as in Section II except that in Section II we did 
not explicitly show the dependence on a forcing function. That is, y = N(x) 
could be defined by y(t) = n(x(t)+f(0) or by y(t) = n(x(t)) + /(*). 

t Actually, an operator related to LN. 
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Cesari" considers the real differential system 

x = g(x, t) 

x - <* , • • • , x n ) (99) 

x, = p,(z, , •■■ ,x n , t) j = 1, •■■ ,n. 

(For the specific conditions imposed on the above functions see Ref. 

14). If 

x t (t) ~ a j0 + X) (a,, cos sat + &,-. sin swi) (100) 

« = i 

and m is a positive integer the vector function Px = (PiXi . . . , P n x„) 

is defined by 

PiXfit) = 0,0 + 2 Co/, cos su« + by.sinaw*) 

i=l, •••,n. (101) 

The operation ff {x) = (Xi ,... ,X n ) is defined by 



Xi(t) = E - (-&.■ cos sw ' + «m sin ■"$ i - 1, • • • 1 n. (102) 

■ am+l Sa > 

The operation F(x) is defined by 

F(x) = H(I - P)g{x), (103) 

where I is the identity map and 

gix) = (g x x, ■•■ , g„x) ^^ 

g f x = gMt), /] 0' - 1, ••• ,n). 

Letting T = P + F, Cesari determines conditions for the existence 
of fixed points of x = Tx. He uses both Banach's fixed point theorem 
(contraction mapping theorem) and Schauder's fixed point theorem 
(which does not give uniqueness but requires weaker conditions). He 
then shows that if y is a fixed point, it satisfies 



If 

then 



h = gMQ,Q + P,<Si- QiV)- (105) 

PMi - 9,V) = 0, j-1, '" ,n (106) 

y(t) = g(,y(i), (107) 
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Cesari discusses the solution of (106) which may be reduced to deter- 
mining a Galerkin approximation. Cesari worked in the space of 
square integrable (periodic) functions and Knobloch 15 adapted his 
approach to the space of continuous functions. 

While the above method has the use of truncations of Fourier 
series in common with our approach, there seems to be a closer rela- 
tionship between Umbo's approach and ours. 

Urabe" considers the real nonlinear periodic system 

where X (a;, t) is periodic in t of period 2tt. If 

xjfi = flu + V2 E a„ cos nt + 6„ cos nt, (108) 

n=l 

a Galerkin approximation* of order m is obtained if one can deter- 
mine the 2m + 1 coefficients a Q} a t ,b x ,. . . , a m , b m that satisfies the 
following equation: 

f -ijr «,.]* 

1 m t *2* 

-\- - y\ i cos nt / X[x m (s), s] cos ns ds 

+ sin nt I X[.v m (s) , s] sin ns ds>- (109) 

"Urabe considers the problem of determining whether there is an exact 
periodic solution near an approximate (Galerkin) solution, .t . He 
determines conditions under which an iteration starting at x converges 
to an exact periodic solution. His proof, while not explicitly mention- 
ing a fixed point theorem, is closely related to the contraction mapping 
fixed point theorem and uses the fact that a contracting iteration 
sequence must stay within a certain sphere centered about the initial 
point.f Our approach is in the same spirit but we take a more general 
viewpoint at the beginning. The basic theorem is derived in an arbi- 



* The method of equivalent linearization is essentially a first-order Galerkin 
approximation. 

tit may be shown that Urabe 's result (Proposition 3, p. 125 of Ref. 16) is 
essentially equivalent to requiring that the operator derivative have norm less 
than one (a contraction) in the appropriate sphere. 
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fcrary Banach space where derivatives of operators are fruitfully 
used. The more general viewpoint is very simple conceptually and 
also permitted the easy extension to difference-differential equations. 
Urabe in Ref. 16, and also in Ref. 17 and 18, considers many aspects of 
Galerkin's method for differential equations not touched on in our 
study. Also, see his comments on Cesari's method on p. 121 of Ref. 16. 

VIII. CONCLUDING REMARKS 

The development of analytical methods (other than asymptotic 
methods) for the equivalent linearization technique with autonomous 
systems remains a very important area for investigation.* Whether a 
modification of the contraction mapping theorem (perhaps using a 
subspace) might be applied to this problem remains to be seen. In con- 
nection with autonomous systems, a question perhaps more important 
than the one we have considered (if equivalent linearization indicates 
a periodic solution, does there actually exist one?) is the following: If 
a non-trivial periodic solution exists, will the method of equivalent 
linearization indicate it? A typical engineering use of the describing 
function is to determine conditions under which no self-sustained oscil- 
lations are predicted. The engineer would like these same conditions to 
also imply that there are no oscillations in the original (exact) system. 
Urabe 10 has shown that the existence of a periodic solution will (un- 
der certain conditions) imply the existence of a Galerkin approxima- 
tion of sufficiently high order. The equivalent linearization technique 
is essentially a first-order Galerkin approximation and the first-order 
approximation may not be high enough to indicate the existence of 
a periodic solution according to a result of the type of Ref. 16. It 
would be very useful to determine conditions that would answer the 
question. This question is related to that rased by Aizerman's con- 
jecture. 

Leaving the problem of autonomous systems we find our adaptation 
of the contraction mapping theorem to be quite convenient in ana- 
lyzing equivalent linearization in forced systems. The calculations are 
straightforward and require no difficult mathematical argument in 
the execution of the basic idea. It is hoped that the method may prove 
useful in justifying and refining approximations. 

It should be clear that our approach is easily adapted to the dual- 



* A theory of autonomous systems, due to Urabe, is outlined in Chap. 3 of 
Halanay. u 



2432 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1967 

input describing function approximation (see, e.g., Gibson, 10 p. 402).* 
The remarks in Section VI concerning difference-differential equations 
apply in that case also. That is, the only essential difficulty is in ob- 
taining the dual-input describing function solution (which has nothing 
to do with our method of investigating the accuracy of such a solu- 
tion). It should be noted that the dual-input describing function 
method has been used primarily for two sinusoids with commensurate 
frequencies (one an integral multiple of the other) and is actually 
equivalent to a Galerkin approximation. When the ratio of the two 
frequencies is irrational, we are in the realm of almost-periodic func- 
tions where analysis can get much more complicated. Boyer has 
presented an interesting approximate method of analysis (an account 
of which is given in Gibson, 19 p. 408ff.) for an input consisting of two 
sinusoids with incommensurate frequencies but with one much larger 
than the other. Analysis of this method would be of interest. 
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APPENDIX A 

Guide to Some Mathematical Background Reading for Engineers 

The purpose of this appendix is to make reading suggestions to en- 
gineers interested in this work but who are not familiar with the 
mathematics used. The primary reference on functional analysis used 
for this work is Kantorovich and Akilov. 6 A more elementary intro- 
duction is given in Kolmogorov and Fomin.°f Lucid introductions to 
the theory of differential equations are given in Hochstadt 13 and 
Struble. 20 The theory of oscillations is extensively covered in Minor- 



* The dual-input describing function was apparently first used by J. C. West, 
J. L. Douce, and R. K. Livesly. In Ref. 7 there is an example of the existence of 
a subharmonic solution to Duffing's equation. This is actually an example of the 
dual-input describing function. 

t The reader should be cautioned that some of the terminology is not stand- 
ardized among American and Russian writers. For example, Kantorovich and 
Akilov do not require a compact set to be closed while most American authors 
do. Also, a linear operator is necessarily bounded according to Kantorovich and 
Akilov but not necessarily bounded according to most American writers (Kolmo- 
gorov and Fomin's definition agrees with American writers on this point). 



EQUIVALENT LINEARIZATION 



2433 



sky 1 which also has a discussion on difference-differential equations. 
Kaplan treats Fourier series and finite Fourier transforms on an ele- 
mentary level. 4 Further information on Fourier series (but still on an 
elementary level) can be found in Tolstov. 21 

APPENDIX B 

Derivatives in a Banach Space 

The following material is abstracted from Chap. XVII of Kantoro- 
vich and Akilov. 6 

Let P map an open subset 12 of a Banach space X into a subset A 
of another Banach space Y. Let x t 9, and suppose that there exists 
a linear* operation U mapping X into Y such that for every x t X 

P(x + te) - P(x ) 



inn 

l->0 



= U(x). 



(HO) 



The linear operation U is said to be the derivative of the operation P 
at the point x . We write this 

U = P'(x ). (Ill) 

The derivative thus defined is the Gateaux or weak derivative and 
U(x) is the Gateaux differential. 

If the convergence relationship of (110) is satisfied uniformly with 
respect to all x e X with || x || = 1, then the operation P is differentiate 
at the point x and the derivative P'(x ) is called the Freohet or strong 
derivative. 

APPENDIX C 

To discuss the satisfaction of (59) , let 



Z = 3 I b I CT 



Consider a, f, and co fixed with / ^ and let a t [0, 1). To show that 

lim Z = 

6-0 



A 


+ 


CT 


| bA 3 |] 
4 


1 


— ex. 



we must show that 



lim I bA 2 I = 



* Kantorovich and Akilov 8 include boundedness in their definition of linear. 
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since 



Z = 3CT\\ bA 2 I +— : + 



~2 I bA <T ■, bA 



I 2 |3l 



1 - a ' VI - a) 16 

From (42) we have that 



hA2 = i{A + »- a 



Also 



so that 



REFERENCES 



/ 



lim A = 



lim bA 2 = 0. 

J.-0 



Several of these works refer to still other works which contain helpful informa- 
tion. For example, see the first footnote on page 2407. 

1. Minorsky, N., Nonlinear Oscillations, D. Van Nostrand Co., Inc., Princeton, 
New Jersey, 1962. . 

2 Truxal, J. G., Automatic Feedback Control System Synthesis, McGraw-Hill 
Book Co., Inc., New York, 1955. 

3. Graham, D. and McRuer, D., Analysis of Nonlinear Control Systems, John 

Wiley and Sons, Inc., New York, 1961. 

4. Kaplan, W., Operational Methods in Linear Systems, Addison-Wesley, Bead- 

ing, Mass., 1962. . 

5. Bass, R. W., Mathematical Legitimacy of Equivalent Linearization by De- 

scribing Functions, in Automatic and Remote Control, ed. J. F. Coales, 
Butterworths, London. 1961. 

6. Kantorovich, L. V. and Akilov, G. P., Functional Analysis in Normed Spaces, 

(translation), Pergamon Press Limited, Oxford, England, 1964. 

7. Holtzman, J. M., The Use of the Contraction Mapping Theorem with Deriva- 

tives in a Banach Space, unpublished work. 

8. KrasnoseFskii, M. A., Topological Methods in the Theory of Nonlinear In- 

tegral Equations, (translation), Pergamon Press Limited, Oxford. England, 

9. Kolmogorov, A. N. and Fomin, S. V., Elements of the Theory of Functions 

and Functional Analysis, Vol. 1, Graylock Press, Rochester, N. Y., 1957 
and Vol. 2, Gravlock Press, Albany, N. Y., 1961. . 

10 Sandberg, I. W., On the Response of Nonlinear Control Systems to Periodic 
Input Signals, B.S.T.J., 43, May, 1964, pp. 911-926. 

11. Halanay, A., Differential Equations: Stability, Oscillations, Time Lags, 

(translation) Academic Press, New York, 1966. _ 

12. Stoker, J. J., Nonlinear Vibrations in Mechanical and Electrical Systems, 

Inte'rscience Publishers, New York, 1950. 

13. Hochstadt, H., Differential Equations, Holt, Rhinehart and Winston, New 

14. Cesari, L., Functional Analysis and Periodic Solutions of Nonlinear Differen- 



EQUIVALENT LINEARIZATION 2435 

tial Equations, Contributions to Differential Equations, Vol. 1, No. 2, 
_ J. Wiley and Sons, Inc., New York, 1963, pp. 149-167. 

15. Knobloch, H. \V\, Remarks on a Paper of L. Cesari on Functional Analysis 

and Nonlinear Differential Equations, Michigan Math. J., 10, 1963, pp. 
417-430. 

16. Urabe, M., Galerkin's Procedure for Nonlinear Periodic Systems, Arch. 

Rational Mech. Anal., 20,1965, pp. 120-152. 

17. Urabe M. and Reiter, A., Numerical Computation of Nonlinear Forced Oscil- 

lations by Galerkin's Procedure, J. Math. Anal. Applic, U t , 1966, pp. 107- 
140. 

18. Urabe, M., Periodic Solutions of Differential Systems, Galerkin's Procedure 

and the Method of Averaging, J. Diff. Equat., 2, 1966, pp. 265-280. 

19. Gibson, J. E., Nonlinear Automatic Control, McGraw-Hill Book Co., Inc., 

New York, 1963. 

20. Struble, R. A., Nonlinear Differential Equations, McGraw-Hill Book Co., 

Inc., New York, 1962. 

21. Tolstov, G. P., Fourier Series, (translation), Prentice-Hall, Inc., Englcwood 

Cliffs, New Jersey, 1962. 



